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Abstract 


A numerical calculation for the full thermal dynamics of colliding nuclei has been 
developed. Preliminary results are reported for the thermal fluid dynamics in such 
processes as Coulomb scattering, fusion, fusion-fission, bulk oscillations, compression 
with heating, and collisions of heated nuclei. 

I . Introduction 

The preceding paper by Ray Nix summarized high-energy nuclear collision calculations 
using relativistic fluid dynamics. In this paper we shall discuss lower-energy nuclear 
collisions, and we shall change the emphasis from predicting outgoing multi-particle 
distributions to attempting to calculate the interior thermal properties of the nuclear 
fluid during and after a collision process. We present this study as an example of a 
methodology for treating the dynamics of any hot compressible fluid drop with known 
material properties. In our particular case we have arrived at the material properties of 
our exotic nuclear fluid by solving a meson field equation integrally with the thermal 
dynamics. For more common materials the thermal fluid dynamic methods presented here can 
be used in a more straightforward manner with an analytical equation of state. 

The original motivation for this study came, not from nuclear physics or from fluid 
dynamics, but from a critical problem in astrophysics: as increasingly more accurate 
modelling is being included in calculations of stellar collapse, it is becoming more 
difficult to get those calculations to result in a supernova explosion (1). The 
calculations end, not with a bang, but with a whimper. Supernova explosions, while not a 
common sight, are definately known to occur, of course, and the astrophysicists amongst us 
would like to be able to predict when they occur. 

A major remaining uncertainty in stellar collapse calculations is the equation of state 
of nuclear matter (NEOS) . The equation of state is a fundamental relationship between the 
pressure and the density and internal energy of a fluid, in order to predict supernova 
explosions, the current collapse calculations of J. Wilson seem to demand rather radical 
departures from the currently fashionable NEOS. 

The material description of a star nearing the end of its collapse is given in Table I. 
At densities around lO^^ gm/cm^ the constituent nuclei are almost touching and the 
temperature of the nuclear material is around 8 - 10 MeV. The heat capacity of the 
nucleus will determine whether the nuclei can thermally disintegrate under these 
conditions. Thermal disintegration is important because a free nucleon gas creates extra 
thermal pressures that help to generate an explosion. 

Heavy nuclei, in the laboratory at least, exhibit a large number of excited states, and 
their extrapolated heat capacity appears to be so large that the nuclear matter is too 
"stiff" during collapse, producing non-exploding stars. While we don't understand all the 
implications at this time, J. Wilson's collapse computations might be interpreted to mean 
that supernova explosions require a "softer" NEOS at high densities and temperatures than 
current nuclear theory predicts. 

Until recently, laboratory heavy -ion experiments were confined to nuclear matter near 
zero temperatures, and extrapolations to higher densities and temperatures led to widely 
varying results. During the last decade, however, high-energy heavy-ion accelerators have 
allowed nuclear researchers to observe experimentally the collisions of heavy nuclei and 
their subsequent excitation and breakup. The goal fror the NEOS point of view is to 
obtain some information about the heating and compress n of the composite nuclear fluid 
before it disintegrates. In heavy ion collisions excitation temperatures as high as 80 
MeV are not uncommon (as deduced from the Maxwellian distribution of the emitted 
fragments; see Figure 1) (2). At low excitation energies the temperature dependence is 
that of a Fermi gas and at high excitations the temperature follows a classical Boltzmann 
gas law. High compressions have not yet been definitively observed, though there is 
intense interest in creating high compression because of a theoretically-predicted phase 
transition at about twice normal density. At this density the nuclear matter might 
exhibit a pion condensation, a state which contains many bound pions. Even higher 
compression might lead to a "density isomer", a stable state of very high density (3). 
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Thus it is extremely important, for the astrophysicist as well as for the nuclear 
physicist, that we understand the NEOS at hiqh densities and temperatures. One way to 
approach this understanding is to adopt a NEOS and study its consequences for heavy ion 
collisions, looking for predicted phenomena that will help to pin down the exact form of 
the NEOS. We have begun such a study. In the next section we describe the NEOS approach 
taken by J. Wilson for our initial studies. 
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Table I; Stellar collapse 
conditions just before bounce. 



Figure 1: Maxwellian distribution of 
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rom the heavy-ion collision 
U at 250 MeV/amu. 


II . The Equation of State 

There are at present two approaches generally used to represent the NEOS: first, the use 
of effective two-nucleon potentials such as the Skyrme interaction (4); and second, a 
direct calculation of the equilibrium energy per nucleon in terms of the underlying meson 
field such as the Walecka NEOS (5). All nuclear fluid dynamics (NFD) calculations up to 
the present have used the first approach; in this paper we report preliminary results from 
an effort to use the second approach in conjunction with NFD. 

The Wilson nuclear equation of state used in this work was derived by solving the 
classical meson field equations using a Lagrangian energy density given by 

^ H 1^) e e ^ 3 ^^ 

where ^ is a scalar field with effective mass ra*“m, (I+ o/Pq) • T represents 
the electrostatic field, and pg and g are adjustable parameters. The density of 
nuclear matter is p, the internal energy per nucleon is c, and Ze is the electric 
charge per nucleon. This choice of Lagrangian was motivated mostly by a desire for 
simplicity since our lengthy two-dimensional calculations include both dynamic and thermal 
effects. Also it was felt chat the astrophysical calculations for which this model was 
intended contained greater uncertainties in other quantities than would be introduced by 
this simple but partially ad hoc Lagrangian. 

Whereas the standard approach begins with a micro.jCopic Lagrangian and eliminates the 
quantum meson fields by replacing them with classical averages, Wilson’s approach begins 
with a combination of microscopic fields and macroscopic densities which is designed to 
reproduce most of the known properties of nuclear matter. In particular, by including the 
density-dependent effective meson m*, he retrieves in a simple manner the saturation 
properties induced by the repulsive vector field in more elaborate models. 

Wilson's approach is related to the Thomas-Fermi class of approximations used in nuclear 
and atomic theory. The use of this unusual means of arriving at the equation of state has 
to be justified by a comparison with more refined NEOS models and with experimental data. 
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Out NEOS parameters are determined by fitting the dimensionlesa variables q, oo> 
m so that they reproduce the binding energy, the measured density profiles, and the period 
of compressible radial oscillations ("breathing mode"), for the complete spectrum of 
nuclear masses. The fitted values are g*0.66, oo"^‘32, and m«0.73. It is a 
significant achievement that a model with these three constants alone can reproduce all of 
the above data. Figure 2 shows the resulting binding energy vs. mass curve compared with 
experimental values. Figure 3 compares the predicted density profiles for b and 
‘^Fe with the experimentally deduced profiles. The equation of state is plotted in 
Figure 4. The compressibility implied by this NEOS is 200 HeV, which is close to the 
experimental value of 210 MeV. For comparison we also show in Figure 4 the 
state>of-the-art (sero-temperature) NEOS of Halecka (S) as well as the sxyrme NEOS used by 
Stocker et.al. in their recent studies of nuclear collisions (6). The latter implies a 
rather high compressibility of 29S MeV, and the Walecka NEOS is even n‘- itter at higher 
densities. The Wilson NEOS used in our studies is relatively soft, which may be more 
appropriate for supernova applications. 



Figure 2; Comparison of experimental 
binding energies with Wilson's 
predictions for all masses A. 



Radkm (fm) 

Figure 3: Density profiles for Fe and 
Pb predicted by the Wilson nuclear 
equation of state. 



Figure 4: Comparison of current 
equations of state. 


The work presented in this paper represents 
ongoing studies of dynamic heavy-ion collision 
processes. These studies compare our predictions 
using the above Wilson NEOS with experimental 
data and with the results of other dynamic 
models. Our main goals are to evaluate the 
consequences of this NEOS and to explore the 
effects of temperature, density and pressure 
gradients on the collision dynamics, both for 
cold and hot nuclear matter. These calculations 
are unique in that they simultaneously treat the 
full macroscopic thermal and dynamic fluid 
processes, including shock waves, while directly 
solving a classical scalar meson field equation 
for the nuclear binding. While the results shown 
here are for central collisions of inviscid 
identical nuclei, we also have the capability of 
treating central collisions of any two nuclei 
with classical viscosity. 
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III. Thermal Physics 

Many effective interactiona used in current NPD calculations include compressional but 
not thermal effects. The compressional part of the internal energy is usually represented 
by a three or four-term expansion of the internal energy as powers of the local density. 

In our calculations we treat both thermal and compressional contributions to the internal 
energy. The compressional contributions are calculated directly from the meson field 
Lagrangian as a function of the internal energy, as described above; the thermal 
contributions are calculated as described below. 

Thermal effects are included in these calculations by directly computing thermal 
pressures, using the thermal energy of a nonrelativistic Fermi gas: 

Here Q is a tensor von-Neumann artificial viscosity. The pressure is obtained from a 
Y-5/3 law: 




This equation is 
by 


valid for all temperatures. The degeneracy energy a 



and 





.ssure are given 


(4) 


(5) 


in units where I ■ 10"13 cm, t ■ 1/3 x 10"23 sec, and m ■ 1 amu (the velocity of 
light is unity) . 


The temperature is not an explicit variable, so we obtain it from the internal ener , / 
through a temperature-defining relation given by: 

^ I + cT^ 




( 6 ) 


The constants a, b, and c were derived 
limit of a cool Fermi gas (x-kT/ejeg) ' 


theoretically to give the correct expression in the 

€ a I J5.TV X 
10 


the correct expression in the limit of a hot Boltzmann gas: 

^4W = I (8) 

and the exact thermal energy in the transition region where the chemical potential is 
equal to kT. 


Thus this temperature-defining relation has been fitted to vary smoothly from low to high 
temperature regimes. The appropriate numerical values are a - 2.96, b - l.S, and 
c ■ 0.66. 
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TV. Computational Method 


To solve the full thermal dynamical problem, J. LeBlanc has employed sophisticated 
hybrid Euler ian-Lagranqian finite difference techniques developed over the years at th? 
Lawrence Livermore National Laboratory. Three conservation equations are solved 
simultaneously in an Eulerian frame: 

(1) the equation of cont uity 

~ -ijfiL 4 I • U, 

dt a-t ^ 

(2) energy conservation 


^ = -ev®- = 


( 9 ) 


(10) 


(3) momentum conservation 


d ^ U*V^U 

dt d“t 




^V(i eZfj 


( 11 ) 


where Q is the von Neumann artificial viscosity for shocks. The fielcs ^ and y are 
solved by iterating two Poisson's equations: 


(4) meson fi 


2 2 

field: ^ ^ - rY7^ where ■+ 


( 12 ) 


(5) electrostatic fi 




(13) 


These five equations represent the full thermal dynamics of this compressible system. 

In our calculations the fundamental variables are the pressures; the internal energy is 
calculated from Pth'^V, while the zone accelerations are calculated from Ptotal* '^hich 
includes the degeneracy pressure. Leaving out jme technical details, we could describe 
this process as t-.<inq y-law excursions off the T»0 degeneracy adiabat. 

Figure 5 displays the sequential operations performed by the computer code. The 
Eulerian mesh typically contains around 2000 nodes (the dots shown in the graphics). The 
basic Lagranwian parts (the pressures and accelerations at time t) are calculated first. 
Next the Eulerian terms are calculated using a second-order accurate implic.t scheme on 
two staggered time-levels. Finally the meson and electrostatic fields are calculated 
using an overrelaxation method. The whole ca'''ulation is extremely fast, taking on the 
average 60 uaec/cycle on a 7600. The time step 6t is calculated according to several 
stability and accuracy algorithms, and then the code cycles on to calculate everything for 
time t 6t. The output is edited in the form of graphics and files on fiche or in 
movies if desired. 

In particular we edit plots of momentum vectors and contours of density, temperature, 
pressure, and the accelerating potential. Examples are displayed in subsequent figures. 
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HEAD-ON COLLISION TRAJECTORIES 



’^’'Ture ‘j; Flow chart for 
advanced Eulerian thermal 
hydrodynamics code of J. LeBlanc 


a) Coulomb tcatunnt: 




(b) Compotita fittion: 


Ic) Fusion window: 
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(d) Fusion-fmion: 
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Figure 6; Possible features of 
symmetric heavy ion collisions as 
a function of projectile energy. 


V. Nuclear Collisions 


In this section we present a variety of collision phenomena predicted by this model, 
along with comparisons with other dynamic nuclear model predictions and with experimental 
data. At the lowest bombarding energies we do not expect our results to be accurate 
because we neglect predominant quantal shell effects. Low-energy collisions are subsonic 
and basically incompressible. When we increase the relative projectile velocity to the 
nuclear sound speed, around v/c = .16, the collisions become supersonic and then 
interesting compressible and thermal effects are predicted. Of course in a detailed 
thermal calculation the compressibility 


K 



(14) 


depends on the local state of the matter, and therefore so does the sound speed 



A£ 


If local conditions are such that Cg ' smaller than in the stiff ground-state nuclear 
matter, then interesting compression'., -re possible. 


( 15 ) 


In our collision calculations we are seeking detailed knowledge about such thermal 
compression effects. Our hope is that with detailed thermal understanding we will have 
more control over predicting the compression and dynamics than has heretofore been 
possible. We shall now review generally what happens in these collisions as the 
projectile energy is increased. 
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ORIGINAL PACiZ 15 

OF POCFi QUALITY 


1. Head-On Collision Trajectories 

Figure 6 illustrates head-on collisions, in the center-of-mass system, of two identical 
nuclei. Five general types of process are possible, in order of increasing energy. If 
the incoming kinetic energy is so low that the long-range Coulomb repulsion absorbs all 
the initial kinetic energy before the shorter-range nuclear forces can fuse the nuclei 
together, then the nuclei turn around and accelerate backward, as shown in trajectory 
(a) . At slightly higher energies the attractive “proximity force" (which is like the 
molecular Van der Waals force) can partially fuse the nuclei. If the initial kinetic 
energy is used up before a stable “saddle shape" is reached, then the repulsive Coulomb 
forces predominate and the partially fused system separates ag.^in, as in trajectory (b). 

At even higher energies it is possible, at least for the lighter nuclei, to fuss the 
composite system inside the saddle shape (where attractive nuclear forces just balance the 
repulsive Coulomb forces), at the same time absorbing most of the initial kinetic energy. 
This results in a stable fused compound nucleus as in (c) . At medium-high energies the 
nuclei fuse to form a flat disk, which then rebounds out to an elongated spheroid. Here 
the kinetic energy remains so high that the system fissions on rebound, nroducing the 
fusion-fission process of trajectory (d) . These higher-energy collisions can result in 
some exotic shapes. Finally, at relativistic energies the nuclei tend to splat into an 
expanding disk whose disintegration was discussed by the previous speaker. 

The existence of a “fusion window" (c) between composite fission (b) and fusion-fission 
(d) depends on the size of the compound system. Theoretically a compound nucleus with no 
kinetic energy reserves will undergo fission if its “fissility" 

X “v 0.02 z2/a = Coulomb energy/surface energy (16) 

exceeds unity. For head-on collisions this occurs roughly for identical projectiles 
heavier than Sm. Thus we would predict that colliding Ca nuclei would have a fusion 
window but U nuclei would not. The fusion window was of great interest a decade ago 
because it made the generation of "superheavy* elements possible. Subsequent experiments 
implied that the fusion window was much narrower than predicted, probably because of 
viscous dissipation during the collision. We plan to study viscous effects, but the 
simulations presented here are for inviscid nuclear matter. 


Figure 7 shows the results of our Coulomb scattering calculation of Ca on Ca at an 
initial code kinetic energy of 0.5 MeV/amu in the CMS. The maximum temperature rema 
relatively cool at 2 MeV and the compression (p/oq where Dq = 0.145 amu/fm3) 
remains near the normal nuclear value. Our computer runs place the fusion barrier 
somewhere between 0.5 and 0.6 MeV/amu, while experimentally it lies at 0.66 MeV/amu. 


ins 


(7) 


Dtraity 


t-142 


284 I 



Momantum 


10 


20fm 0 



f-Figure 7: 


40 


Ca 


+ 40 


E„„=0.5 

cm 


Ca 

Mejj 

amu 



20 fm 



Momantum 






274 



0R5G1MAL AA23 s?5 
OF PCCR Q'JALir/ 


Dmiity 




0 10 20 fm 


Monwntum 




Dtraity 


Momentum 



40 40 

Figure 9: ’"Ca + ”Ca at E =1.5 MeV 

cm 

amu 


Figure 10: 


40 


Ca + 


40 


Ca at E =5 MeV 
cm ' _ 
amu 


That is, we currently underestimate the barrier, perhaps indicating a little too much 
nuclear binding. Time-dependent Hartree-Fack (TDHF) calculations, on th? other hand, tend 
to overestimate the barrier to between 0.7 and 0.9 MeV/amu, depending on the choice of 
iteraction force (7). 

A Coulomb scattering calculation for a 239y collision is shown in Figure 8. Notice 
the interesting sequence of distorted shapes during closest approach. 

Figure 9 shows a fusion process (Ca on Ca) at 1.5 MeV/amu. The composite system necks 
into a prolate spheroid and then flattens into an oblate disk, from which it rebounds into 
a spheroid without a neck. This fused compound nucleus then oscillates very much like a 
classical charged liquid drop, as shown by the later frames in Figure 9. During the 
initial shock of fusion the compression reached a maximum of 1.3 and the average 
temperature was around 5 MeV, indicating strong local heating. The temperature tends to 
increase quickly during the initial fusion process and then to build up slowly thereafter 
(see Figure 14) . 

Next in this sequence. Figure 10 shows a high-energy Ca on Ca collision at 5 MeV/amu 
which demonstrates fusion followed by fission on the rebound. Here the oblate disk gets 
very extended but the nuclear forces are able to pull it back so that it shoots out into 
an extended prolate shape which soon fissions. The maximum compression at impact was 
1.4. During the expanded oblate stage the average compression was sub-nuclear at 0.7, and 
the average temperature cooled down to 3 MeV. 

If this collision is simulated at even higher energies, then the oblate disk expands so 
much that it does not stop, and the composite system is presumed to disintegrate into 
small particles. 

Thus our Ca on Ca calculations display the appropriate behavior as a function of initial 
kinetic energy. The dynamic fusion window for Ca on Ca appears to occur for Ecm between 
0.6 MeV/amu and 4 MeV/amu for the case of this low-viscosity nuclear fluid. The lower 
limit agrees with the expei^imental of 0.66 MeV; the upper limit has not been measured to 
our knowledge. We have not yet established our upper limit precisely; it is somewhere 
between 3 and 5 MeV/amu. 
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ORlGli^AL PACZ IS 
OF POOR QUALITY 

Figure 11 displays the collision and subsequent fission of a 239u system at EQg^sl.S 
MeV/amu. The accelerating potential contours are plotted here. The collision results in 
a ternary fission of the composite. At lower energies a binary rission results from 'he 
collision. 
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A rinal example of collision dynamics is shown in Figure 12. Here ISOsm collides with 
itself at 5.2 MeV/amu in the center of mass. The interesting feature is the torodiai 
shape that forms from the oblate disk. Similar toroidal shapes have been observed in TDHF 
calculations. We also see them in higher-energy 239u collisions. The torus expands out 
with low compression and then collapses in on itself with high compression. The composite 
then fissions on the rebound. 

2. Composite Nucleus Oscillations 

During and after a low-energy fusion collision, the composite nucleus exhibits two kinds 
of oscillation. During the collision it undergoes a series of ra(.id energy oscillations 
associated with compression waves in the heated composite. After the fusion is complete 
the heated compound nucleus oscillates in a roughly classical liquid drop manner. These 
processes are exhibited in Figure 13, where we show the total kinetic energy as a function 
of time for the collision of two Kr nuclei at 0.5 MeV/amu in the CMS. 

The period of the surface oscillations after fusion is 18 x 10"22 sec, which is close 
to the period predicted by the classical expression for frequency of surface oscillations 
of a charged liquid drop. Thus at these lower energies the thermal effects generated 
during the collision do not greatly affect the overall classical kinetic behavior of the 
compound nucleus. This predicted oscillation behavior can also be taken as evidence that 
our code is working reasonably at low energies. 

At higher collision energies a stable compound nucleus is not formed because it fissions 
on rebound. The heated fission fragments, however, also undergo roughly classical 
oscillations as they move apart from each other. 
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Above nuclear excitations of the order of SO MeV, shell effecrs disappear in heavy 
ni.cl»i. This corresponds to a temperature of about 6 MeV using the usual conversion for a 
Fermx gas: (8) 

(Ec.m.-VVnucleon 
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At such temperatures, macroscopic features 
like shapes and fission barriers are influenced 
by the excitation, and heating effects are 
important in medium-and high-energy nuclear 
physics. In collapsing stars the compressed 
nuclear fluid is at temperatures of 2 - 10 MeV, 
so thermal changes in the nuclei comprising the 
material could affect the bulk properties of the 
star. Figure 14 compares experimentally deduced 
fragment temperatures with the theoretical Fermi 
and Boltzmann limits (8) . 
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Figure 14: Nuclear Temperatures. 


Compression during collisions is of great interest because of the possibility of phase 
changes. At twice the ground state nuclear density a pion condensation has been 
predicted. Figure 15 compares the maximum compressions obtained during a barely subsonic 
^^Ca fusion-fission collision at 5 MeV/amu and a supersonic collision at 19 Mev/amu. 

The outstanding feature is the rapid expansion into an underdense torus for the supersonic 
case. This is accompanied by cooling as shown by the corresponding temperature curves in 
Figure 16. When the torus collapses, a compression of 1.6 results. It is possible that 
careful engineering of nuclear dynamics could produce compressions near 2.0. The subsonic 
collisions remain mostly incompressible, with some compression followed by expansion 
during the fusion process. 
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Thermal expansion also accompanies the higher-energy collisions. At> the nuclei fuse the 
composite heats up and expands. These thermal properties translate into a 
temperature-dependence for such macroscopic guantities as the surface tension, the Coulomb 
energy, and the volume energy, as pointed out by Sauer et. al. (9). The static thermal 
(T»2 MeV) Hartree-Fock calculations of the latter produced smaller changes in the density 
profile than we see in our heated ^®Ca nuclei. 

At high temperatures the fission barrier is lowered due to a decrease in surface 
tension. We plan to study this process in more detail in the future. Conversely, nuclear 
material at high temperatures, as in stellar material, might be expected to exhibit 
different fusion windows and cross sections, both because of thermal expansion and because 
of softer surface density profiles. We will now discuss this interesting possibility. 

4. Collisions of Heated Nuclei ; 

Our results for collisions of hot nuclei are very preliminary.. This type of calculation 
was accomplished by allowing two cold "laboratory" nuclei to execute a fusion-fission 
reaction, and then to let the heated fission fragments impact their mirror images by means 
of reflection boundary conditions on the right-hand mesh boundary. We did this in order 
to achieve realistic velocity, temperature, and potential profiles in the "initial" heated 
nucleus . 


The reasons for studying collisions of hot nuclei are twofold; first, to examine nuclear 
collisions in the astrophysical environment; and second, to explore possible ways of 
creating abnormal nuclear composites using secondary collision fragments in the 
laboratory. The latter should really be studied using heated projectiles and cold 
targets, and we have begun such studies as well. 


A successful fusion of the hot fission fragments from a 239u ternary fusion-fission 
reaction is shown in Figure 17. Unfortunately our code did not edit the masses of the 
independent fragments, but this hot fused composite is evidently very massive, with a 
radius around 8.S fm. This intriguing result should await a careful examination of this 
new type of calculation before it is to be believed. 
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Figure 17: Collision of hot secondary 
fragment with its mirror image (on HHS) . 


The results at** indeed very 
preliminary, but there is a suggestion 
here that the fusion window might be 
enlarged by pre-conditioning the 
projectile to higher temperatures or an 
expanded radius. (However, we have not 
included viscous effects yet, and the 
classical viscosity, which goes to 
t 1/2, will tend to close the fusion 
window.) Surface oscillations of the 
secondary fragment could also be used to 
change the fusion window if the period is 
engineered, for example, to present an 
oblate disk shape to the next target 
nucleus. To observe these phenomena we 
should concentrate on secondary 
reactions, tailoring the primary energy 
to produce the desired characteristics in 
the secondary. 

A phenomenon which may be related has 
recently been reported by Friedlander et. 
al. (10). They observe that some 
secondary fragments display an 
anomalously short mean free path for 
their charge. We would like to suggest 
that these "anomalons" may be secondaries 
from central collisions that were heated 
and then expanded and cooled, resulting 
in a larger geometric cross section than 
would be calculated for a normal nucleus 
of that mass. However, for this 
explanation to be possible the heated 
nuclei would have to stay in an expanded 
state for very long times (10'^ sec). 

This might occur if they cannot radiate 
away their internal energy. We plan to 
investigate this possibility quantitatively. 
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VI. Conclusions 


We have developed a coniputational method for treating the full thermal dynamics of 
nonrelativlstic compressible nuclear matter. In particular we have begun a systematic 
study of central collisions of heavy nuclei. Preliminary results were presented for such 
processes as Coulomb scattering, fusion, fusion-fission, bulk oscillations, compression, 
heating, and collisions of heated nuclei. 

Our main goal at this point was to exercise the calculation through the above variety of 
dif f icult-to-model phenomena to verify that it gives generally reasonable results. The 
predicted behavior in the above processes agrees gualitatively with experimental 
experiences e.q. low-mass projectiles fuse into a stable compound nucleus; high mass 
projectiles undergo a fusion-fission process; low-energy projectiles turn around under the 
influence of Coulomb forces; and perturbed nuclear drops oscillate in a classical manner. 
In producing such overall good behavior we have shown that the Wilson nuclear equation of 
state gives reasonable dynamical results. 

We have barelv begun to analyse the current predictions quantitatively. Much 
information is contained in these calculations. Besides reproducing some known laboratory 
phenomena, we have also predicted some interesting features which may or may not be real: 
temporary toroidal shapes are common in the higher-energy collision calculations, and the 
high-mass composites tend to fission in a ternary manner with a light remnant in the 
center of mass. These ternary processes can be checked experimentally if they persist 
when our calculations are made more quantitative. We also see some very preliminary 
evidence that pre-heated secondary nuclei can fuse more easily than cold primary nuclei. 
This would have important implications in both nuclear physics and astrophysics. 

Our next step will be to try to systematically fine-tune the code predictions of the 
available experimental data to a point where the calculations ate quantitatively 
accurate. Our immediate goal is to calculate the detailed temperature, pressure, and 
density gradients (including shock waves) during both cold and heated nuclear collisions. 
If this can be done, then the opportunities presented to us are enormous. In a long- 
range view, one can envision an era of nuclear fluid engineering in which one can 
carefully control the dynamics of nuclear collisions to produce abnormal composites, 
perhaps with applications totally unheard of today. 
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